o 

(N 



Many-particle correlations in non-equilibrium Luttinger liquid 

I.V. Protopopoyi^^^ p g Gutman^, and A. D. Mirlin^'^'^ 

^ Institut fiir Nanotechnologie, Karlsruhe InsUtute of Technology, 76021 Karlsruhe, Germany 

^ L.D. Landau Institute for Theoretical Physics, Kosygin str. 2, 119334 Moscow, Russia 

^Department of Physics, Bar Ban University, Ramat Gan 52900, Israel 

'^Institut fiir Theorie der kondensierten Materie and DFG Center for Functional Nanostructures, 

Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 

^Petersburg Nuclear Physics Institute, 188300 St. Petersburg, Russia 

We develop an operator-based approach to the problem of Luttinger liquid conductor in a non- 
equilibrium stationary state. We show that the coherent-state many-body fermionic density matrix 
as well as all fermionic correlation functions out of equilibrium are given by one-dimensional func- 
tional determinants of the Fredholm type. Thus, the model constitutes a remarkable example of a 
many-body problem where all the correlation functions can be evaluated exactly. On the basis of 
the general formalism we investigate four-point correlation functions of the fermions coming out of 
the Luttinger liquid wire. Obtained correlations in the fermionic distribution functions represent 
C^ ' the combined effect of interaction and non-equilibrium conditions. 

(N ■ 

I. INTRODUCTION 

1^ ■ It was understood long agc^ that conventional description of interacting fermionic systems in the framework of 

c/3 Landau Fermi-liquid theory is not applicable in one dimension (ID) because of infrared divergences. In particular, 
f~^ the second order of perturbation theory yields a singularity in the self-energy at the Fermi surface, SftS ~ elne, 
implying a vanishing quasi-particle weight. This violates the one-to-one correspondence (which plays the central role 
in the construction of the Landau theory) between unperturbed electronic states and elementary excitations of the 
interacting system. The lack of correspondence is a hallmark of the emerging strongly correlated electronic state - 
Luttinger liquid (LL). 

The physics of LL is known to be relevant for many systems available in experiment. The applications of this 

concept include carbon nanotubes^Ti^, semiconducting and metallic nanowires^, edge states of the samples in the 

Q ■ quantum Hall regime-"^ and spin laddersiSiii 

O . In view of singular infrared behavior, it is hard to access LL by conventional methods of many-body fermion 

perturbation theory. To overcome this difficulty, a powerful approach-the bosonization-was developed. It is based 

on the fact that in ID fermionic creation and annihilation operators have simple representations in terms of bosonic 

fields describing such observable as charge and spin density. Thus, any ID fermionic system is equivalent (as long as 

considered energies are not too high) to some (generally interacting) bosonic system. This Fermi-Bose equivalence was 

^O first discovered on the level of correspondence between correlation functions in fermionic and bosonic theories^—. 

'•ft Later on, underlying operator relations were derived by Haldane^®, providing the solid basis for bosonization (for a 

^j , recent detailed exposition, seai^). 

t — ' Bosonization has proven to be a very efficient tool for tackling ID interacting fermions. In some cases it allows to 

^D \ obtain exact solutions of the problems that are highly non-trivial in the fermionic language. The canonical example 

] ' of this is the Tomonaga-Luttinger model equivalent to free bosons. Even in the case when bosonization does not 

produce free bosons, it often constitutes a convenient starting point for the development of the theory. Using this 

approach, effects of backscattering22i2i, impurities^^, and underlying periodic potential^^ on the LL were explored. 

Many particular realizations of the LL states were addressed, both theoretically and experimentally^i^i. 

Nowadays, there is a growing interest in non-equilibrium phenomena in the LL phase. In particular, in a recent 
C^ ' experiment^, the tunneling spectroscopy of a biased LL was carried out. A similar approach was implemented 
to study experimentally carbon nanotubes^^ and quantum Hall edges^^— . The experimental advances motivated 
theoretical interest to quantum wires out of equilibrium'^'^^'^^ . In a related line of research, non-equilibrium chiral 
ID systems have been studied theoreticallyfiSri*^^ mainly in application to experiments on quantum Hall edge-states 
intcrferometry. 

In a series of papers by two of the authors and Gefen^i^ the theory of non-equilibrium LL was developed. The 
formalism developed in these works combined the bosonization with the non-equilibrium Keldysh action approach. 
It was shown that the calculation of single-particle Green function reduces to the evaluation of certain Fredholm 
determinants, analogous to those encountered in the context of counting statistics^ and non-equilibrium orthogonality 
catastrophe^!^. This allowed to obtain comprehensive results for observables related to the fermionic single-particle 
Green function, including the tunneling density of states, electron distribution function, and Aharonov-Bohm signal. 
More recently, the counting statistics of the charge transfer in a non-equilibrium LL was studied ini^^. 
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However, a number of fundamental questions have remained open. In particular, it is important to understand 
what kind of correlations (if any) are induced by the interaction between the electrons coming out of an interacting 
wire. In particular, are the left- and right-movers that have passed through the wire correlated? We will see below 
that at equilibrium no such correlations exist. On the contrary, in a non-equilibrium LL electrons experience a specific 
type of relaxation. We will show that, in course of this relaxation process, non-trivial correlations in the electronic 
distributions are built. To characterize these correlations is the main task of this work. 

Our approach to the problem is complementary to the one of Ref. 31. Using the operator formalism, we determine 
explicitly the many-body density matrix p of the LL conductor in a stationary non-equilibrium state. When written 
in the bosonic coherent-state basis, this non-equilibrium density matrix p has a form of a Fredholm determinant. 
The density matrix carries full information about many-body correlations in the non-equilibrium LL state. We 
find that the calculation of fermionic correlation functions with density matrix p is greatly simplified by using a 
"refermionization" procedure. In this way, we can evaluate all fermionic correlation functions, with results expressed 
in terms of Fredholm determinants. Having developed our general formalism, we employ it to study four-point 
correlation functions of fermionic fields out of equilibrium and explore the correlations in the occupation numbers of 
the outgoing electrons induced by the interaction. 

The structure of the paper is as follows. Section [H] contains a description of a model of a LL conductor out of 
equilibrium. In Sec. lIIII we give a short review of bosonization and refermionization and introduce concepts needed in 
the main part of the paper. In Sec. IIVI we obtain the many-body density matrix of a LL. Further, we show how the 
bosonization- refermionization allows one to obtain correlation functions of a non-equilibrium LL. We verify that the 
results for the single-particle Green function agrees with that obtained previously in Ref. [Sll In Sec.|V]we explore the 
four-point correlation functions of fermions that emerged out of the LL wire. Our result reveals non-trivial correlations 
in fermionic occupation numbers. Section [VII contains a summary of our findings. 

II. THE MODEL 

Let us specify our model of the non-equilibrium LL conductor. We consider a ID wire populated by spinless electrons 
of two chiralities (labeled hy rj ^ R, L) interacting via local density-density interaction. The wire is connected to 
non-interacting electrodes. To model the later we will assume that the electron-electron interaction is switched off 
outside the central part of the wire (Fig. [T]) so that the Hamiltonian of the problem reads 

Hee = Ho + \j dxg{x) {pI{x) + pI{x)) . (1) 

Here Hq represents the kinetic energy part of the Hamiltonian; Pr(l){x) is the density of right- (left-) moving electrons 
at point X. The function g{x) approaches a constant value deep inside the interval |x| < 1/2 (region II of Fig. [Ij and 
is zero for |x| > 1/2 (regions I and HI). This way of representing the non-interacting leads was used extensively in the 
literature for the analysis of the transport properties of Luttinger liquids^"—. 

The non-equilibrium state of the wire is induced by the injection of the electrons from the leads. Following 
Refs. l3ll we assume that the right- moving electrons coming into the wire from the left lead have a stationary non- 
equilibrium distribution function ^/^(e), while the left- moving electrons coming from the right lead are characterized 
by the distribution function nL{e). The simplest non-equilibrium state arises when the leads are held at different 
temperatures T/j, T^ and different chemical potentials fin, p^ so that the distribution functions of the incoming 
electrons are of the Fermi-Dirac form 

""(^^ - l + e(^-.„)/T, • (2) 

Such a state was named the "partial equilibrium" in Ref.—. More complicated distribution functions n,,, e.g. double- 
step distributions, can be generated if one assumes that the leads are diffusive conductors. We refer the reader to 
the Refi^i for the comprehensive discussion of the possible experimental realizations of the non-equilibrium Luttinger 
liquids. 

Throughout the paper we assume the absence of the electron backscattering in the wire which can be caused e.g. 
by impurities. In order to avoid the electronic backscattering also from the boundaries of the interaction region we 
suppose that the function g{x) is smooth on the length scale of the Fermi wavelength. Withing this limitation, the 
formalism developed in this work can be applied to the case of arbitrary dependence of the interaction constant on 
X. However, the most interesting situation arises when boundaries of the interaction region are sharp on the scale 
of the typical wavelength It* = Vp/T* of the bosonic excitations. Here Vp is the Fermi velocity and T* is some 
characteristic width of the distribution functions n^j (e.g. T* ~ max{Tfi,TL) in the case of the partial equilibrium). 
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FIG. 1: (Color online). Schematic view a LL conductor driven out of equilibrium by the injection of non-equilibrium elec- 
trons (with distribution functions n™(e)) from non-interacting leads. The leads are modeled via the assumption that the LL 
interaction constant g{x) is space dependent and vanishes outside the central part of the wire (|a;| < 1/2, region II). 

Under such circumstance, there is significant scattering of the bosonic excitations on the boundaries of the interaction 
region which is the primary source of the effects we are going to discuss in this paper. For this reason we concentrate 
bellow on this "sharp boundary" limit and model the x-dependence of the interaction constant by 

g{x)=geil/2-\x\) . (3) 

In Ref;^ the single particle correlation functions of the model described above were investigate in great detail. 
While within the interacting part of the wire those correlation functions are highly non-trivial due to the combined 
effect of interaction and non-equilibrium, outside the interaction region the single particle Green functions are much 
simpler, reflecting the free dynamics of the electrons. Specifically, the single particle Green functions of the right 
electrons to the right from the interacting region (region III) has a form of the Green functions of free fermions with 
some energy distribution n^"*(e). In the same way, the left-movers in the region I can be viewed as free particles 
with energy distribution n|^"*(e). The correlation functions of incoming electrons, i.e. right-movers in the region I 
and left-movers in the region III, are not affected by the interaction. The distribution functions n°"* differ from n^ 
(except for the case of complete thermal equilibrium) , signaling the redistribution of the electrons over energies upon 
the passage through the interacting region //. 

Now an important question arises. Is this redistribution (which, crudely speaking, can be referred to as relaxation, 
although the resulting electronic distribution functions are not of the Fermi-Dirac type) is the only effect of interaction 
on the outgoing fermions? In particular, are the left- and right-movers coming out of the wire independent? Note 
that energy relaxation in our device is of a rather special character. Due to the strong constraints imposed by the 
integrability of the LL model, no energy relaxation can occur in a uniform LL. The relaxation actually takes place at 
the boundaries of the interacting wire and involves the scattering of the bosonic excitations from one chiral brunch 
to another. Since bosons are crudely speaking electron-hole pairs, we may anticipate that the relaxation should be 
accompanied by the buildup of correlations in the distribution functions of left- and right-electrons. We show below 
that this is indeed the case and the density matrix of outgoing fermions shows strong correlations. In particular, the 
irreducible correlation functions of the form ((n^(a;, e)n^/(x', e'))) between the occupation numbers of the outgoing 
fermions are non-zero. These emerging correlations are in the focus of the present work. We stress that the correlations 
in question are specifically non-equilibrium effect absent under the equilibrium conditions. 

We are now ready to present the formalism that we use in this work to explore correlations in a non-equilibrium 
LL wire. We start with a short review of the standard bosonization approach (in the form of "constructive 
bosonization" ^^i^^ ) in application to equilibrium interacting fermions. We make particular emphasis on construc- 
tions that will be subsequently employed in the analysis of the non-equilibrium state. 

III. BOSONIZATION AND REFERMIONIZATION 

Let us consider the ID fermionic system with the Hamiltonian 

H = Jdx^P+{x)^^^^^ix) + V,^t. (4) 

Here kp is the Fermi momentum and Vint represents the four-fermion interaction. Within the Luttinger liquid model 
one linearizes the spectrum of the fermions near the Fermi points and adds unphysical (but irrelevant at energies 



smaller than the Fermi energy) states below the bottom of the Fermi sea. This way of treating the interacting 
fermions is justified for the description of low-energy properties of the system. At high energies effects related to the 
curvature of the spectrum (not included in the LL model) may become sizeable, see, in particular Ref. [SOi 

With the bottom of the Fermi sea pushed down to minus infinity, the Hilbert space of the problem is spanned by 
the fermionic creation and annihilation operators aZ'j, , a,j_fc labeled by momentum — oo < k < +00 (counted from 
kp) and the chirality rj = R, L (in formulas, we will also use the notation ?/ = 1 for right- movers and 77 — —1 for 
left-movers). It is convenient to make the fermionic momenta discrete, assuming that the system is placed on a ring 
of large circumference L. The thermodynamic limit L — >■ 00 is taken at the end. The physical fermionic field is 
decomposed into the sum of left- and right-moving parts according to (A — > cxi stands for the high energy cutoff) 

i,{x) ^ Vfl(a:)e''="" + Mx)e~"'^' , (5) 

V',(x) = ^^a,,,e^'=--I'=l/^ (6) 



The Hamiltonian of the LL model acquires the form 

H^VrY^vk {a+^a^,k - r^°(fc)) + ^ ^ dxg{x){pL{x) + pu{x)f . (7) 

h.r] 

Here riJJ(fc) — {—r]k) are occupation numbers of the left- and right-movers in the free-fermion ground state; the 
interaction term was taken in the form of local density-density interaction with space-dependent interaction constant 
as discussed in the previous section. 

The central role in the Luttinger-liquid model is played by the Fourier components of fermionic densties 

Pv,q = Yl K,k+q°"i^k (8) 

k 

which can be used to construct the set of bosonic creation and annihilation operators according to (775 > 0) : 

(9) 





K,q^K',q' ^^m'^qq'- (10) 

Let us denote by \Nii,Nl) the state of the system which is a filled Fermi sea with N^tt^^ extra particles in the 

right (left) brunch, i.e the state characterized by the distribution functions n ?, = Q {—rik + 2ttNjj/L). All these 
states are annihilated by bjj^q and are vacuum states from the point of view of the bosons. Any other state of the 
fermions can be generated by the action of the bosonic raising operators onto \Nji, N^) . Thus, the operators b^^q , &+g 
together with the particle number operators A^^ and the Klein factors Fj^ , i^+ changing the total number of fermions 
of corresponding chirality form the complete operator set. In particular, the free-fermion Hamiltonian Hq can be 
reexpressed in terms of bosons as 

H, = VFj2^{m)\q\b+gb,.q + ^ ^ A^„(^r, + 1) , (11) 

r),q V 

while the fermionic field operators are given by the famous bosonization identity 

V'+(^) = /Je-'^''(-)F+. (12) 

The phase (f'qix) is related to the correspondmg density Pr;(x) via 

Pri{x) = —d-^^nix) . (13) 

ZTT 

Explicitly, in terms of bosonic creation and annihilation operators the density and phase fields read as follows, 



<^,(x) ^^J2Q{m)J^^ {^-""X,, - e^'^^^^'O + ^N,x, (14) 
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p,{x) = )^ e im) \j 2^ (e-'^^<,, + e*^^6,,,) + -N, . (15) 
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FIG. 2: (Color online), a), b) Bosonic wave functions of the in- scattering states outside the interacting part of the wire, c), 
d) Bosonic wave functions of the out- scattering states. 

In Eqs. (|11[ [Til [T5|) the last terms involving the total number of fermions are usually referred to as the "zero mode" 
contributions to the energy, bosonic phase and the density. We keep them in the equations above in order to maintain 
the completeness of our short presentation of the bosonization technique. On the other hand, this contributions do 
not play any significant role in the non-equilibrium effects we are going to discuss below. For this reason, we omit the 
zero-mode terms altogether in the subsequent formulas. 

The power of the bosonization approach rests upon the fact that the interacting fermionic Hamiltonian ([T]) can be 
reexpressed as a quadratic function of bosonic operators. Thus, H can be diagonalized by linear transformation of 
bosons. Bellow we will need explicit expressions for the the bosons diagonalizing H. To derive them, it is convenient 
to start with the equations of motion for the density operators pri{x,t), 



d_fpnix,t)\ d 

dt \PL{x,t) r dx 




PR{x,t) 
PLix,t) 



= 0. (16) 



We decompose the solution of the equation above into the normal modes according to 

Pni 



,(a;,i) = ^e(r;'g)yjr(y;, l^l(^,:,)5+^^e^l^l^-'-K v,|,|(ry,x)5,,,e-^l^l^-*) . (17) 

q,v' 

Here b and &+ are the new bosonic operators, and Utj' ,\q'\ ('7; q) are the coefficient functions. It is convenient to organize 
the latter into two-component columns Uj-i^q{x) ~ (u,,_g(a;, i?), w,,^g(a;, L)) satisfying the equation 

d 

- lVFqUj^,q{x) + — [Tzh{x)Ur,,q{x)] = , (18) 



where t, is the third Pauli matrix and 



Vf 



ajx) g{x) 



M.)^('^- ,J'^J. (19) 

Note that in order to construct the densty fields we need only the solutions of Eq. (fT8|) with q > 0. Thus, in the 
discussion of the properties of Urj.q{x) we will assume that q > 0. On the other hand, for q < one can define w,,^^(a;) 
via the complex conjugate 

Uv,q<0ix) = K^-qi^)- (20) 

This choice is consistent with Eq. (jlSp . 

Imposing different boundary conditions, one can construct two sets of solutions of the Eq. ([T5)) . We denote them 
by {q > 0): u™ and u^'i*. Far from the interacting region the solution ^'"fj.f^.') consits of a right (left) incoming 
wave and the waves scattered of the interacting region. These are the in-states in the standard terminology of the 
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scattering theory. On the other hand u°^ and u°^l represent the out- scattering states. They contain the outgoing 

wave and the waves incident on the barrier. The asymptotic behavior of of u,"g " is summarized on Fig. [21 We 
have introduced the bosonic g-dependent scattering coefficients of the interaction region. Due to the time reversal 
symmetry the transmission amphtude tq for the right wave w™^ coincides with that of the left wave, whereas reflection 
amplitudes for the right and left waves are related via f^/r^ = —V/t. If we assume the boundaries of the interaction 
region to be sharp on the scale of the relevant bosonic wavelength and adopt the model of the rectangular-shaped 
interaction constant ([3]) we flnd: 

^1 = 1 _ 7^2g2.if,Z ' (21) 

^-i{l-K)ql 

TR^q = rL,q^-2tn^— ^^j-^-^ Sin Kql. (22) 

Here K is the LL parameter of the interaction region, while TZ = {1 — K)/{1 + K) and T = \/l — TZ"^ = 2\fK jiK + 1) 
are reflection and transmission amplitudes of a single interaction boundary. 

Clearly, outgoing solutions w°"* are not independent of u™ . The two sets of solutions are related by the elements 
of the bosonic scattering matrix of the interaction region Vq and tq. 

'^R,q ~ ^q^R,q ~^ ''^L,q'^L,q ' (23) 

<i = ^k.<:. +*;<:.• (24) 

Choosing u™ (a;) as the coefficient functions in Eq. (J17p . we obtain the decomposition of the density in terms of 
the in- bosons 6™ and &+™, 



q'.ri' 

Comparing (|25l) to ([T5|) . one can read off the relation of the initial and in-bosons: 

KT = E ©W) [uM,,W,v'K,q, - n\^a~q',v')b,r,q] . (26) 

ri',q' 

Operators 6™ are defined via hermitian conjugation and uV^„{q') are the Fourier components of u™ (x): 

<,r,{'i'.v') = \\rj\ f dxe-^^'^ul-^ix,r,') . (27) 

By construction, operators 6™ and &+™ diagonalize the Hamiltonian H: 

if = T/^^e(r;g)M6™+6™ ■ (28) 

The connection (J13l) of the phase fields (p^ to the densities can now be used to derive the relation between (p^ and 
the incoming bosons. 



^''(^) = 4 E ®(^''^'h/ :^ {K'^wi^^^^Hw - <"k'i('7:^)&;",') • (29) 

q' ,ri' V ' ' 

Here we have introduced the two-component columns vj^q{x) given by (the matrix h{x) was defined in (fTO)) ) 

V:;^q,ix)^h{x)u';^^q,ix). (30) 

Since tp^{x) ^ exp [— i(^^(a;)], equation (|29p connects the fcrmionic flcld to the in-bosons. 



For the purposes of this work, it will be convenient to perform refermionization of the in- bosons 5™ by introducing 
new fermionic fields "01" (a;) and their Fourier components a*\ according to 

V'+^lx) - 4f E <re"*'" = -^e-^'"(-)F+ , (31) 



VI V " VI^ 




^l"{x) = z^ e(7yg)^/— (e-''-5+;" - e^^^b" ) . (32) 



The inverse relation reads 



Since the in- bosons solve the LL Hamiltonian, the in- fermions a*"^. are just free fermions with the time evolution 

«;':feW = «;':fee-'"=^-*. (34) 

Comparing Eqs. ([^ and (15^ and taking into account the boundary conditions imposed on u™^, we see that for 
rjx — ?> — oo (e.g., for right fermions at a; — > — cx)) the physical fermionic field ip,i{x) is identical with "iptix) , 

ipTj{x) — ip^{x) , 77X — >■ — oo . (35) 

So far, we were focussing on the in- solutions u™ . The out- solutions w°"* allow us to define the out- bosons 5°"* 
and fermions a°^^ by equations analogous to (^5]) . (PT|) . and (15^ up to a replacement m — ;> oui everywhere. These 
out- fermions and bosons also solve the LL Hamiltonian. Just as the solutions u™ and u°"* , the bosons 6™ and 5°"* 
are related by the bosonic transmission and reflection coefficients. 



^L,q — ^R,q"R.-q + ''q'^L.q 



°L.q - ^R,q"R.-a^'^q"L.a ■ ('^ ' > 



The out- fermions represent the physical fermions after the crossing of the interaction region in the sense that (cf. 
Eq.dSSl)): 

0„(a:)=07*(x), r;a;^+oo. (38) 

The refermionization procedure described above turns out to be the crucial step in the solution of the non-equilibrium 
LL model. The in- fernfions introduced in pip will provide us with the starting point for the construction of the 
non-equilibrium density matrix of a LL wire. The out- fermions will be useful in our discussion of the correlations 
among the electrons leaving the wire. 

IV. NON-EQUILIBRIUM LUTTINGER LIQUID 

A. Density matrix 

In the previous section we have summarized the bosonization representation of the LL Hamiltonian. Let us now 
turn to the analysis of the non-equilibrium LL in the setup shown on Fig. [T] First, we need to understand what is the 
(many-body) density matrix p describing the situation of Fig. [TJ This density matrix should satisfy two conditions: 

a) it should be stationary with respect to the LL Hamiltonian (O; 

b) any correlation function of the fermionic fields (-0^^ (xi, ti) . . .^^\xi,ti) . . .) evaluated with the density matrix p 
should be identical to the same correlation function of free fermions with the density matrix 



Po = -^ cxp 



rj.k 



(39) 
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as soon as all coordinates of the fermionic fields satisfy ijiXi < —1/2. In Eq. (l39l) Z stands for the normalization factor 
and the parameters e^ (fc) are determined by the distribution functions of the fermions coming into the wire via 



n^(fc) 



1 



1 + e''i(''') 
Both conditions a) and b) are fuUfiled by the density matrix 



(40) 



1 



ri,k 



(41) 



The condition a) is satisfied due to the equation of motion ([34| for the incoming fermions. The condition b) is satisfied 
due to the fact that the fermionic fields iprnixi) are just identical to Tpir'{xi) as soon as rjiXi — >■ — oo. 
In equilibrium, e^(fc) = rjkVp /T , the density matrix in terms of bosons is simply given by 



Peq ^ exp 



' T 



Y.^{m)\<iK7K, 



(42) 



Together with Eq. (|29p this implies that the evaluation of the fermionic correlation functions in equilibrium is reduced 
to the evaluation of the averages of the type (griiQ) and g^{q) are c- valued functions) 



Z[9ri,9v] = (exp 



^'d{m)9v{Q)f} 



in+ 



exp 



Y.^{m)g;{q)K, 



(43) 



with the gaussian weight (|42t , which is a straightforward task. 

Translating the general non-equilibrium density matrix ()4ip into the bosonic language is a much more non-trivial 
task. As shown in Appendix |Al it turns out to be possible to evaluate the matrix elements of p in the basis of bosonic 
coherent state |iV, /?) (eigenstates of the bosonic annihilation operators 6™ ) in the form of one-dimensional Fredholm 
determinants. Here we state only the final result of this analysis, referring the reader to Appendix VK\ for the precise 
definitions and details of the calculation: 



Dn = det 



{P*^,PlNR,NL\p\NR.NL.PR,fiL) ^ DR{P},,fiR)DL{PhPL) . 

1 - e*''(^)n^''(fc)e-*''(^)(l - n„(fc)) - e-*^(=^)(l - n^''(fc))e*^(=^)n^(fc) 



(44) 
(45) 



N„ 



In the determinants ()45p the distribution functions n,,(fc) and riri^ {k) are considered as operators diagonal in the 
momentum space. On the contrary, operators $j)(x) are diagonal with respect to the conjugate variable x and given 

by 



$„(fci-fc2) = e(rKfci-fc2)). 



27r 



L\ki 



7Pr],ki—k2 



(46) 



The density matrix (j44|) , (|45|) , along with the Hamiltonian ([28l) and the expression for the fermionic operators ((29|) 
contains the whole information about the problem in the language of non-interacting bosons. The natural next step is 
to calculate n-point fermionic correlation functions (that determine various physical observables) . This will be done 
in Sec. HVBl 



Electronic correlation function 



Our goal now is to evaluate many-point fermionic Green functions that have the form (|43p , with the bosonic density 
matrix given by Eqs. (j44 |) . (j45 |) . While this can be done by a direct calculation in the bosonic language, such a way 
turns out to be quite tedious. A shorter way is to perform a refermionization of the bosons 6™. Indeed, according to 

. Further, the density matrix is quadratic in 



m+ 



Eq. ((33)) . 6*"+, &™ are quadratic functions of the in- fermions a 

terms of the in-fermions as well, see Eq. (PT|) Thus, the average (1^5]) can be expressed as trace of a certain operator. 



which is an exponential of a quadratic form with respect to fermions a*" . The evaluation of the trace leads to (see 
Appendix [B] for details) : 



Z[g*,g]^An[6R{x)]AL[6L{x)] , (47) 



Ajj[6jj{x)] = det 



(l - n°(fc) + e-^^-^^^nO/fc)) ' (l - n,,(fc) + e-'*"(")n^(A:)) 



(48) 



6,{x) = ^ ^ J^Mvq) (.9.,,e'^^ - 9;,,e~^n • (49) 
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m 



In Eq. (j48p we can recognize a one-dimensional functional determinant of the type discovered in"^^ in the context 
of a single-particle Green function. The first factor in square brackets of Eq. ()48p involves the zero-temperature 
distribution function n^.{k) and serves as a regularization. It ensures that at zero temperature Z[g*,g], which is an 
average of a bosonic normal-ordered expression, is identically equal to unity. 

Having derived the general result, we turn to evaluation of many-point fermionic correlation functions. Let us 
consider the average of a product of n operators ■0,,. {xi,ti) ,i — 1 . . .n and n operators V-'it {xi, i,) , i = n + 1 . . . 2n, 

*^»)l...r;2„(2;i,^l, • ■■X2n,t2n) = Uprni^litl) ■ • • V'r)^ (a;™ , ^n) V'^„+i (a;„+l , i„+l) . ■ ■ i^^^^{x2n-,t2n)) ■ (50) 

Representing the fermionic fields by bosonic exponets according to ([T2l [29| and brining the product of the exponents 
into the normal ordered form we get 

-^V---')2„(2;i,^l,---a;2n,i2n) = M^^^^,^^^(xi,ti,...X2n,t2n)Z[g*,g] , (51) 



Here (^i = 1 for i = 1, . . .n, while Ci = — 1 for i = n + 1, . . . 2n. In the expression ([FT]) the first factor Af" arises 
due to normal ordering of the bosons and is nothing but the corresponding correlation function of ^-operators in 
zero-temperatures LL. This factor depends in particular on the odering of '0-operators in (|50p . On the other hand, 
the second factor in (j51l) accumulates the effect of the distribution functions of the incoming electrons n^(e). This 
factor equals unity at zero temperature and is the same for all the fermionic correlators which differ only by the 
ordering of the Fermi fields. Explicit expression for the factor Z[g*gr^ is given by the Fredholm determinants pS)) 
with the phases 5,-i{x) expressed in terms of the functions Vri^q{x,ri') via [cf. Eqs. (|49I52I) ] 



^J.i-)-'-^J^±OvZiv. 



^.)g-^gV^t,-»^?x^ (53) 



Equations ([51]) . (|Tf| and ([55)1 together with (IT51[5(H) provide the full solution of the non-equilibrium LL problem. It 
is straightforward to check that the two-point (single-particle) correlation functions found previously in Ref.'^^ are 
correctly reproduced. 

Expression for the phases 5.q{x) acquires particularly simple form in the case when all the fermionic fields in (j50p 
are taken outside the interacting part of the wire, i.e. \xi\ > 1/2 for alH = 1, . . . 2n. Taking into account the boundary 
conditions imposed on the functions u™ we get 

WJ^{x) ^27rf]f^Yl Oe^^f"'"'-''") [6 (r,,x,) (t,<5^,^, + r^,,<5^,-^J + e(-r,,a;,)<5^„J ■ (54) 

Here we have introduced the light-cone coordinates Ui = Xi — rjiVpti. According to ([SI]) , in the case of the electronic 
correlations at the input of the central part of the wire (i.e. for all rjiXi < ~l/2) the phase Sri{x) are completely 
independent from the properties of the interaction region. This ensures the coincidence of the correlation functions 
with that of the free fermions with the density matrix (P^ . On the other hand, for the correlations of the fermions 
at the output of the interacting region (i.e. rjiXi > 1/2) the phases (5,, (a;) are determined by the bosonic transmition 
and reflection amplitudes. 
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V. CORRELATIONS IN THE OUTGOING FERMIONS 

A. Single-particle Green functions 

We are now in a position to turn to the correlations between the electrons going out of the interacting part of wire. 
The simplest quantities characterizing this correlations are the four-point correlation functions of the type (j50l) where 
the right fermions are "measured" to the right of the interaction region, while the left electrons are "measured" to the 
left, i.e. we consider the behavior of Eq. (|50p for ijiXi > 1/2. Intuitively, the electrons which have left the interaction 
region are just free electrons. We can make this statement mathematically precise by recalling Eq. (pS)) showing that 
for rjiXi > 1/2 physical electronic fields coincide with the fields of out- fermions, having just free dynamics 

^^(x,t)=V7*(x,t)=V7*(x-r/Ffi), vx>l/'^- (55) 

Thus, the question arises, if any of the correlation functions (j50l) are non-trivial. However, apart from the dynamics, 
there is another important ingredient of the correlation functions. This is the density matrix. 

At equilibrium, the density matrix in terms of the in-bosons is given by (|42|). Using the interrelation (|37)) of in- 
and out- bosons and refermionizing the out-bosons, one readily finds 



1 
Peg = -^ exp 



rp 7 y 



^^F(<r<fc-"°w) 



(56) 



Thus, at equilibrium the out-electrons have both free dynamics and trivial density matrix. In this situation none of 
the correlation functions (j50|) with T]iXi > 1/2 bears any trace of the interaction in the central part of the wire. In 
particular, the prefactor Af " entering Eq. (j5ip for the correlation functions at the output of the wire is just the free 
fermion zero temperature correlation function. 

The situation changes drastically in a non-equilibrium system. The density matrix p becomes now a complicated 
function of the out- fermions, incorporating non-trivial many-particle correlations. It is important to distinguish 
between the correlation functions which are determined solely by the dynamics and the correlations also governed by 
the electronic distribution. The retarded and the Keldysh single particle Green functions are the simplest examples of 
the former and the latter. For the non-interacting electrons G™(a;, r) ~ 8{x — r/VpT). Thus, for the retarded function 
we immediately get 

G^'^i^x^r) = Gf{x,T)^n\bR = ^\^l\5l - 0] = Gf(x,T) . (57) 

On the other hand, the Keldysh component of the Green function (say for right electrons) G^ixR., tb/2\ x^, —tb,/2) = 
—i[ipri{xR,Tj^,/2),ijj^{xR, —Tfj/2)] = G^ {tj^} receives non-trivial contribution from the combined effect of interaction 
and non-equilibrium, 

G^ixR,Tn/2,XR,^TR/2) = G'^°{0,tr)Ar[6r]Al[6l] , (58) 

where the phases (5,, are given by 

S/J„ix) = -innj f ^e*9(-H-'7-) sin ^^ {tgS.,,R + tl.^S^.l) . (59) 

J iTT I 

Our choice for the arguments of the fermionic fields in the definition of the Keldysh Green function is somewhat over 
complicated, since G^ is actually independent of xr ( as long as xr > 1/2). However we keep xr explicit for the 
convenience of the forthcoming disscussion of the four-point correlations. 

Exploiting the explicit expression ([2^ for the transmition and reflection amplitudes in the "sharp boundary" model 
we get the phases ^,, in the characteristic form of sequence of rectangular pulses of length vpTpr^. The amplitudes of 
the pulses are shown on Fig. [3l Their positions are given by 

XRm ^XR-l + (2m + l)Kl , m = , 1 , . . . (60) 

for 5r{x) and 

XLm — ~xr, + / — 2viKl , m = , 1 , . . . (61) 



for the left phase. Analytically 



5,^{x) ^2TT'^aj^^rnWrR{x- Xn^m)- (62) 
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5ij/27r 



r^ 







r^v? 






r^n^ 





xm 



xri 



Sl/27t 



XR2 



n 



Xl2 



Xl\ 



2-773 



r'Ti 



Xlo 



T'n 



FIG. 3: The phases 5r{x) and 5l{x) governing the Keldysh Green function of the right-movers. The coordinates a;,,.o ,a;,,,i 
are given by Eqs. (|60I61[) . 



Here the coefRcients 



and we have introduced 



2yr}2m 



OLKm ^T TZ 



OiL' 



n, m = 

_^2^2m-i rn^ 1,2, 



Wt{x) ^ Q \ -— " \x\ ] signr. 



(63) 



(64) 



Equation (1581) can be interpreted in terms of the fermionic occupation numbers at the output of the wire, 

1 - 2<,"*(e) = tVp I dTRC^'^'^G^iTR) . (65) 

Since the n°"*(e) are different from the distribution functions of the incoming electrons n,j(e) (except for the equihb- 
rium) , electrons experience relaxation upon crossing the interaction regiori^. 

B. Two-particle Green functions 



Let us now turn to a deeper characterization of the density matrix for the outgoing electrons that is provided by 
the four-point correlation functions. Particularly intriguing are interaction-induced correlations between the left- and 
right-movers. To reveal them we consider— 

Mrl{xb^,tr;xl,tl) ^ { ^b^[xr,y) ^"^Ry^R^ — ^) V'l (a^i, y j , V'l (^^l, - yj ) + Gf (rflJGf (tl) . (66) 

were xr > 1/2 and xl < —1/2; the square brackets [•, •] stand for the commutator. Note that (1551) is the only non- 
trivial irreducible correlation function of two right and two left fields. Replacement of one of the commutators in (|66p 
by anticommutator would immediately lead to the decoupling of left and right fermions under the average. It is easy 
to see that the Fourier transform of (j66p has the meaning of the irreducible correlator of two fermionic distribution 
functions, 



{{nR{xR,€R)nL{xL,eL)}} = ^ I dTRdTLM{xR, tr; xl, TL)e 



„i<^RTB.+i<^LrL 



(67) 



Each of the phases (5,,(x) corresponding to (1661) is a sum as of a pulse sequence encountered previously in the 
discussion of the Keldysh function G^ and an analogous contribution originating from the pair of left operators in 
(j66| . This last contribution consists of pulses of the length Vftl located at 



''Rm 



-XL- 1 + 2mKl , 



m 



0,1 



(68) 
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for Sii{x) and 

x'L^^XL+l-{2m + l)Kl, TO = 0,1,... (69) 

for the left phase. 

The dependence of (IS71) on the coordinates Xrj and the "optical length" of the wire Kl can be easily understood. We 
note that the scale lT*/yF for the times r^ supporting non-zero correlation functions is set by the inverse characteristic 
width of the energy distributions of the incoming electrons. Under the assumption that the length of the interacting 
wire is large and Kl exceeds It* we can neglect the mutual effect of the non-overlapping pulses in ^B.iL) ^nd represent 
the corresponding determinants A^[(5,j(a:;)] by a product of determinants for individual pulses. One immediately 
concludes that occupation numbers of left- and right- movers are uncorrclatcd unless the pulses coming from right 
and left operators overlap. The overlap happens if Xj^x = ^r + ^l is close to odd multiple of Kl. We can now recast 
(p7| into the form 

{{nR{xR,eR)nL{xL,eL))) ^ ^ f^''{eR,eL, Xrl - mKl) . (70) 

iriGOdd 

The functions /^^(ej^, e^, y) arc independent of I and decay on the distance It* ■ They are given by 
f,n^{eR, €L,y) = -^ / dTLdTLG^{TR)G^{TL) (ArAl - 1) leRTR + Ieltl , 



(71) 



Mrr{xr,tr;x'j^,t'j{) = ( ijR {xr,y) 'V4 {^R' — §) i'R (^^'"^) ' ^^ (^«' — 2 



^ ^ TT Afl [aR.n-rhWrnix) + aL.nWrL jx + v)] ,^^. 

-'■-'; AR[aR^n-rhWrj,{x)]AR[aL,nWrr,{x)] 

n— max(0,rn) . j u ; 

A _ TT Al [aR,n-rhWri^{x - y) + aL,nWrji{x)] 

-'■-f; /\L[aR,n-rhWrL{x)\AL[aL,nWrii{x)] 

n— inax(0,m) 

Here rh = {m + l)/2, and we use the notations introduced in Eqs. (|551[M1) . 

Similarly to (1661 [67]), one can define the irreducible correlation function of the occupation numbers of the right- 
movers: 

{{riRixR, eR)nRix'j,, e'j,))) = ^ / drRdr'^^MRRixR, tr; x'j,, r;,)e'^«^«+"'«^« , (74) 

+ G^irR)G^Xr'j,). (75) 

In order to be able to interpret (j74p as the correlator of the distribution functions, we assume below that the two pairs 
of -0-operators in ([75)1 are inserted far from each other, so that Xrr = xr — x'^ » It* ■ Under this assumption, just 
as in the case of R-L correlations, (1751) is the only one non-trivial correlation function of four right electronic fields. 
Consideration of phases S,^{x) corresponding to ((75)) shows that the ([71)) is non-zero provided that Xrr « 2Klm and 
can be decomposed as (we set m = m/2): 

{{nR{xR,eR)nR{x'j,,e'j,)))= ^ f^^''ieR,eL,XRR - mKl) , (76) 

meeven ,m^o 

ineR, eL,y)^-^J dTRdT'„Gl{TR)Gf{TL) {ArA^ - 1) e---«+-«^« , (77) 

_ -^ AR[aR^n-'mWrR{x) + aR.nWr-^{x + y)] 

-'■-'; AR[aR,n-niWrR{x)]AR[aR^nWr'{x)\ 

^ _ jT Al [aL,n-,nWrAx) + aL,nWr'^{x - y)] 

^^ Al [aL,„-rnWri (a;)] A^ [aLnWr'ix)] 

n— inax{0,?n) x. v / j l ' n ^ ' j 

The term with to = is omitted from the summation due to our assumption that Xrr > It* ■ Again, investigation 
of the functions fm^i^R, £l,x) requires evaluation of the functional determinants A^. 
The physics of obtained correlations is discussed in Sec. IV CI 
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C. Discussion 

We assert that the correlations discovered in Sec. IV Bl represent a quantum interference effect. To clarify this point, 
we make a digression and consider a very simple case of one right boson with momentum q populating the wire. This 
situation is described by the density matrix 

Pib = b]^-\0){0\b^^^, (80) 

where |0) is the ground state. Using the connection of the in- and out-bosons and then performing the refermionization 
we translate pib into the out-fermions 



2tt 
Tq 



P^b-'-T. (^.4r+,«?el + rR.^AT-.'^T}) |0) (0| {t>tT-,<k + r*R,,al%^ar]) . (81) 



fe.fe' 



We can now evaluate the correlation function of the left and right distributions directly in the fermionic language and 
get 



{{nR{xR,eii)nL{xL,eL))) = -<5"-fl(efl)5"-L(eL) 

+ ^ «+ - nl_) «+ - <_) (4,,i,e''(^-+^-) + r«,,re-^^-) . (82) 

Here lo = qVp is the frequency of the boson in the system; n"^ = 0(— (e^ ± ci;/2)) and 

Snnien) = nTieR) - n°(efl) = —\t,\' {2n\ei,) - n^eR + u) - n"{eR - c.)) , (83) 

Lq 

SnU^L) = nT\eL) - n^en) = ^\rR,g\' (2n"(ei) - n°(6i + uj) - n°(6i - uj)) . (84) 

l^q 

(85) 

The two terms in (|82)) have distinct physical origin. The first one originates from the probabilistic nature of the 
boson transmission-reflection process. It would remain unchanged upon the replacement of the density matrix by the 
statistical mixture 

~Pib = \t,\'b+';^^0){0\bZ + \rRj'b+^^l\0){0\brt, . (86) 

This term favors anti- correlations of electrons at energies of equal signs. This is a typical anti-correlation between 
the mutually exclusive events. Indeed, suppose the boson we have injected into the system was transmitted through 
the interacting part of the wire. In this case at the output of our device there is a particle-hole excitation in the right 
branch, i.e. a right electron at positive energy and a right hole at negative energy. At the same time the left branch 
has no excitation. On the contrary, if the boson was reflected there is a left electron at positive energy and a left 
hole. The right branch is empty. Thus, in such a probabilistic description there is no way to have an excited electron 
both in the left and right branches, which results in anti-correlations between the fermionic occupation numbers at 
energies of the same sign. 

Note that the probabilistic contribution to the correlator does not depend on the coordinates where the the occu- 
pation numbers are measured. This is quite general. Suppose we have some complicated density matrix p in terms of 
in-bosons. We translate it into the out-bosons. It is easy to see that if we now switch to the classical description of 
the system by dropping all the non-diagonal elements of the density matrix (cf. transition from ([8T|) to ([86])) we will 
end up with the correlator of the occupation numbers independent on Xri^ . Only the matrix elements of p involving 
momentum transfer from left to right branch can provide such a dependence. 

The second contribution to the correlator (15^ is due to the fact that on the quantum level the scattering of a 
boson creates coherent superposition of the state with a boson in the right branch and a state with a boson in the 
left branch. Thus, the exited particle-hole pair is (virtually) simultaneously present in both branches. To elucidate 
the effect of this quantum term, let us consider a slightly more general density matrix which is a statistical mixture 
of dSOD: 

9 q 
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Now the quantum contribution is modified accordingly, 

// / ■, f ^^^ ,^2^. V^ 27r . g w Q \ sin Kql sin qXB.L 
{{nR{xR, tR)nL{xL,tL))) quantum = "47 7^ 2^ c,— [ur^ - nR_) [n^^ - n^_ j — o 
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'L?^ "+ "-^ ^ "+ "-^ \l^^l2^2^KqU 



We have used here the exphcit expressions (j22p for the transmission and reflection amplitudes in the sharp-boundary 
model. Let the coefficients Cq be peaked at some q = qq. Under the assumption that the peak width is much larger 
than 1/Kl and the energies £r(l) are not too close to qoVp/'i, we can average the expression under the sum over fast 
oscillations oscillations on the scale 1/Kl. The result is non- vanishing only if Xrj^ is close to odd multiple of Kl, 
Xrl — mKl + X, m (£ odd, in which case 

47rT^7?.™ sign m ^-^ Stt , „ n \ / n n \ 

{{nR[XR,eR)nL{xL,eL)))quanturn = ^ ^2 ^'^'^JT ^'^^+ ~ ^R-) \^L+ -n^_)cosqx. (89) 

q 

This result is in agreement with the general coordinate dependence (I7D|) of the correlator of fermionic distributions. 
The density matrices pn, and Pmix are very simple as they contain just one fermionic excitation. The truly non- 
equilibrium density matrix (j4ip . when written in terms of bosons, is much more complicated. It contains infinitely 
many terms representing multiple-boson processes. To collect properly their contributions into the correlation function 
of left and right occupation numbers is the task accomplished by the functional determinants A^[J^(a;)]. In general, 
both classical noise and quantum interference contribute to the result. The lesson we learnt from the analysis of 
simple density matrices allows us to identify clear manifestations of quantum effects. First, this is the dependence of 
the correlation functions on the coordinates, see Eq. (|70l) . Second, these are positive correlations between distribution 
functions of left and right movers at energies of the same sign. We will see in Sec. IV Dl below that the correlation 
functions of a non-equilibrium LL do show both these features. 

D. Partial equilibrium 

A detailed analysis of the functions /^^ describing the correlations in electronic distributions requires a careful 
investigation of the determinants A,J^^(a:)] with given distributions of the incoming electrons. In general, these 
determinants can not be evaluated analytically, except for asymptotic long-time behavior'^^ (which is not sufficient 
for our purposes here). Thus, one has to resort to numerics in order to achieve the comprehensive understanding of 
the correlations in the fermionic occupation numbers. 

Relegating such a numerical analysis for a future work, we focus below on the case of partial equilibrium that can 
be treated analytically. This is the situation when distribution functions of the incoming fermions are of the Fermi- 
Dirac form but with different temperatures Tr and Tj, of the left- and right- movers^^. The functional determinants 
A,j[S,j{x)] become under such circumstance exponentially quadratic functionals of phases^, 

lnA„ = -i- r^qiB,iqVF)-l)\S,{q)\\ (90) 

where Bn{uj) — cotha;/2T,, are the bosonic distribution functions. 

One can now use ([50)1 to evaluate the products of the determinants (|72l73l78l79p analytically. Alternatively, we 
can apply Eq. (|90p to the full functional determinants governing the correlations in the electronic distributions. The 
Fourier components of the phases 5^ are easily read off from (|59p . For example, in the case of of the correlator of left 
and right distribution functions we have 



|c , ^|2 (4^)' 



r 



\5L{qr = 



g2 



\tq\ sm — V\rq\ sm — V2itqr^sn\ — - — sm — - — s\nq{XR+XL) 

,2 • 2 qypTR I ,2.2 qypTL ^,^ „* „.„ qVpTR ._^ qVpTL 



\tq\ sin — hlr^l sin — 2zigr* sin — - — sin — - — s\i\q{xR + XL) 



(91) 
(92) 



Plugging this into (|90l) , averaging the expression under the integral over fast oscillations of the bosonic transmition 
and reflection amplitudes (which is equivalent to the neglect of the interference of non-overlapping pulses in 5,-i{x) 
discussed in the previous section) and using the standard equality 



+ 00 



dx 



-(l-cosax)(coth--lj=ln-j— , (93) 
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FIG. 4: (Color online). The dependence of the correlator {{nR{xR,0)nL{xL,0))) on the sum of coordinates xr + xl- Solid line 
corresponds to Tr/Tl = 10 while the dashed line is the result for Tr/Tl ~ 2. Energy dependence of the correlator of left and 
right distribution functions at xr + xl = Kl (blue dot) and xr + xl = ~Kl (brown dot) is shown on Fig. [5] 
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FIG. 5: (Color online). Energy dependence of functions fi'^{eR,€L,0) and f^i{eR,eL,0) (plots (a) and (b) respectively). All 
energies are measured in units of the largest temperature Tr = IOTl ■ 



one finds that the correlation function of left and right occupation numbers has indeed the form (|70p with 



fm {^R,<^L,x) 



-^ I dTiidTLC 



leLTL+ieRTa/^K 



G^irR)Gt{rL) 



1 _ ch TTTB.(rn + TL) -i _ ch ttTl (tj;-tl 



ch 2TTTitx 



ch2TTTz,x 



1 _ chTrTLJTR+TL) 1 

ch2-!TTLX 



c\i2-kTrx 



(94) 



Here the exponents 7™ are given by 7™ = sign77i7^7?.l™l/(l + TiP'), and the Keldysh Green functions of the outgoing 
fermions are^ 



G'^iTR) ^ 
Gf{TL) 



1 



T^VpTR 

1 

TlVpTL 



^TrTr 1 


■^ V r 

l + TL^ 


shTiTRjR 

\ nTLTL ' 


sh ttTltl 



t^Tltr 



sh t:TlTr_ 
sh ttTrTl 



(95) 

(96) 

(97) 



The Fourier transformation in (j94p can be evaluated numerically. The results are exemplified on Fig. |4l[5]and|6l To 
produce the graphs we have chosen the LL parameter K = 0.2 such that the transmission and reflection amplitudes 
T, TZ are approximately equal. Figure H] shows the zero energy correlator {{nR{xR,0)nL{xL,0))) in its dependence 
on Xr + Xl- Solid and dashed lines correspond to strong {Tr/Tl = 10) and comparatively weak (Tr/Tl = 2) non- 
equilibrium. One observes the characteristic peaks at xr + xl — (/2m + l)Kl. Occupation numbers of left and right 
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FIG. 6: (Color online). Correlator of the distribution functions of left and right electrons integrated over the xr + xl (see Eq. 
I). The ratio of temperatures was taken to be Tr/Tl = 10. 



movers are anti-correlated for xr + xl > and correlated at xr + xl < (of course, the situation will be reversed if 
one assumes that Tr <Tl). 

Figure ([5]) demonstrates the functions f^f{€R,eL,0)- The correlations reach maximum near eR = e^ = 0. In the 
case of fi'^ the anti-correlations are somewhat more extended in the direction eR — — e^ while the correlations in f^^ 
prefer to develop for eR, and e^ of equal sign. This fact leads to a remarkable structure in the "mean correlations" of 
the distribution functions of left- and right- movers as given by 



h{€R,€L) = y- j d{. 



Xr + XL){{nR{xR, eR)nL{xL,eL))) 



(98) 



(The prefactor was introduced to make h dimensionless.) The function h{eR,eL) is shown on Fig. [5] We observe 
that on average the left- and right- movers are correlated when their energies have equal signs and anti-correlated 
in the opposite case. Note that this result is opposite to the one that would be except on the basis of the classical 
consideration of the bosonic transmission-reflection process, see Sec. IV CI 

A similar analysis can be performed for the case of the correlations in the occupation numbers of right electrons 
along. One finds the correlation function of right distributions in the form of (j76p with 






T/2 



fl'- 



j dTRdT'j,e''^^-+''''^''G^{TR)G 



^(r'R) 



-*- ~ ch27rTi,a; -*- ~ ch2TrTnx 

ch7rTn(rH+rjj) ch7rTi,(TR-i 

^ ch 2-nTRX ^ ch2TpTZx 



k) 



- 1 



The exponents 7^ are given by 



7™ = r'7^'™l/(l + 7^2) 



(99) 



(100) 



The analogous correlator of the left distributions can be obtained via simple exchange Tr o Tl. Provided that 
Tr, > Tl all the functions /^^ turn out to be positive, while the corresponding functions f^^ are negative (cf. Fig. 
[7] showing /2^^(eii;, e^, 0) and /|'^(eL, e'^,0)). Thus, the hotter electrons have positive correlations built into their 
distribution, while the colder ones are anti-correlated. 



VI. CONCLUSIONS 



In this paper we have developed an operator approach to the non-equilibrium LL. Using bosonization and refermion- 
ization techniques, we have explicitly determined the many-body density matrix of the system and derived the 
fermionic correlation functions in terms of Fredholm determinants A^. Let us note that usually the correlation 
functions of the many-body interacting system can only be found approximately by truncating the Bogoliubov-Born- 
Green-Kirkwood-Yvon chain. The model considered in this work constitutes a remarkable example of a many-body 
problem where all the correlation functions can be evaluated exactly. 
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FIG. 7: (Color online). The energy dependence of the correlators {{nji{eji,xji)nR{e'fi,x'fi))) and {{nL{eL,XL)nL{^'L,x'j^))) for 
Xn~x'n = 2KI. While the occupation numbers of the hotter (right in this case) electrons are correlated, the occupation numbers 
of colder electrons are anti-correlated. 



We have employed our technique to study the four-point correlation functions of the electrons coming out of the LL 
wire. While the dynamics of the outgoing electrons is free, the corresponding density matrix is highly complicated. It 
incorporates the correlations ((n^(a;^, e^)n^'(a;^', t^'))) in the electronic distribution functions caused by the scattering 
of the LL bosons at the boundaries between the LL wire and the non- interacting leads. The spatial dependence of these 
correlations can be deduced from the general analysis of the corresponding functional determinants A^. It displays 
characteristic interference picks (or dips) of the width It- at distance Kl (the "optical length" of the interacting wire) 
one from another, see Fig. 2] 

For the case of partial equilibrium we have evaluated the Fredholm determinants governing the correlation function 
{{njj{xrj, eri)nrj> {xri' , eri'))) analytically. We have found a non-trivial spatial dependence of correlations, both in the 
occupation numbers of electrons of the same chirality and between the occupation numbers of left- and right-movers. 
Within the "hotter" chiral branch, the positive correlation are developed in the electronic distribution. On the other 
hand, anti-correlations are seen in the distribution of colder electrons (see Fig. [T]). The sign of the correlation 
function {{nji{xii, eii)nL(xL, £l))) depends on the sign of xr + xl- On average, the occupation numbers of right- 
and left-movers are correlated at energies of the same sign and anti-correlated at energies of opposite signs, see Fig. 
ini The obtained results indicate that quantum-interference effects contribute crucially to the correlation functions. 
Intrinsically quantum correlations between parts of the system which can not be understood in a classical framework 
are referred to as entanglement^. In recent years it was recognized that quantifying the entanglement between 
subsystems of a many-body quantum system can provide a clue to many relevant properties of the system, see e.g. 
Ref. 54.55|. Such quantification for a mixed state of a quantum system is in general a highly complicated task^^. It 
remains to be seen to what extent the correlations found in this work are relevant in the quantum information context. 

We stress once again that the correlations studied are a genuine non-equilibrium effect absent in equilibrium LL. 
These correlations are an experimentally relevant quantity and can be measured in a specifically designed tunneling 
experiment. For example, to access the correlator of left and right distribution functions, one can imagine the 
following setup. Suppose that two tunneling probes are attached to the wire at points xn, xl satisfying rjXri > 1/2. 
We can also assume for simplicity that the right probe allows only the tunneling of right electrons, while only left 
electrons can tunnel through the left one. Then (under the assumption that the tunneling densities states in the 
probes are not flat) the correlator of right and left tunneling currents {{lR{V]i)lL{VL))) should be sensitive to the 
correlations in the distribution functions. The voltages T/^ biasing the probes will control the corresponding energies in 
{{nfi{xji, eR)nL{xL,eL)})- Clearly, a more sophisticated analysis is needed to extract the information on the correlator 
of the distribution functions from the correlation in the tunneling currents in a realistic setup. Our analysis shows, 
however, that the correlator of the distribution functions is in principle a measurable quantity. 

Concluding the paper, we briefly discuss possible extentions of the present work. First, under the general non- 
equilibrium conditions, a numerical analysis of the functional determinants A,, is needed to fully understand the 
correlations in the electronic occupation numbers. Second, the formalism developed here can be extended to cover 
some non-stationary states of the LL. For example, in view of the recent experimental developments^, it is very 
interesting to investigate a LL wire exposed to on-demand coherent single-electron sources. In particular, in the 
context quantum information processing, especially intriguing is the entanglement generated by interaction between 
two electrons injected by single-electron sources from the left and right leads. Third, upon a proper modification, our 
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technique should also be usefull for the investigation of the fractional quantum Hall edge states out of equilibrium. 
In this systems it is very interesting to look for manifestation of the fractionally-charged quasiparticles in the non- 
equilibrium correlation functions. 
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Appendix A: Density matrix of non-equilibrium fermions in the bosonic representation 

In this appendix, we transform the density matrix (j4ip into bosonic representation. Throughout this Appendix we 
will be dealing with the in-fermions a™ and in-bosons &™ only. In order to simplify notations, we omit the index "in" 
for all the operators. Further, we focus on right-moving in-fermions. The contribution of right-movers is obtained 
analogously; the total density matrix is the product of contributions of left and right in-particles. 

We are thus looking for the bosonic representation of the statistical operator 



P = ^ cxp 



fe 



eik){a+ak-n\k)) 



(Al) 



Here Z is the normalization factor, n^{k) = 0(— fc) is the ground-state distribution function and e(fc) determines the 
distribution of the electrons via 



n{k) 



1 



1 + e<^('=) 



(A2) 



A particularly convenient basis in the Hilbert space of chiral fermions is provided by the bosonic coherent states 
which are the eigenstates of the bosonic annihilation operators. Each coherent state is labeled by the total number of 
fermions N and the set of eigenvalues /3q of the operators bq. Explicitly 



|iV,/3)=exp 



lq>0 



\N) 



(A3) 



The coherent states (|A3p form the overcomplete basis with the resolution of identity given by 



1-E 



n dl3;d(3q 



N " U>0 

The overlap of two coherent states is given by 



exp 



q>Q 



iiv,/3)(iv,r 



(A4) 



(M,/3*|iV,/3)= V^exp 



El/3.1 

.q>0 



(A5) 



We are interested in the matrix elements of the statistical operator p in the basis of the coherent states. They are 
given by (obviously, p is diagonal with respect to the total number of fermions) 



{P\N\p\N,P) = {N\e^p 



Y.f^:^'^ 



.q>0 



pexp 



lq>0 



|iV) 



(A6) 



We remind that the bosons are proportional to the Fourier components of the fermionic density: 



K 



|^E4+,«.' ^. = y|^E4-,«fc. 



(A7) 



Thus, the operator under the average in (IA6|) is exponentially quadratic in fermions. 
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It is convenient to introduce new fermions: 



Ck^ak{l-n''{k))+a+n^{k), 



„+ — „+ 



= a+(l-n^(A:))+afcn^(A:) 



(A8) 
(A9) 



Here n^{k) — 8 (— fc + %-^) is the distribution function of a Fermi sea with N extra particles. The idea behind the 
transformation (|A9|) is that the state |A^) is nuhified by ah the operators Ck, which simphfies the derivation. In terms 
of the new fermions c we have 



q>0 fci ,fc2 fci ,fc2 

The matrices Uki,k2 s-nd Vk^^k2 earring two momentum indices are given by 



\T^N 



U = $(l-n''')-$-'n 



(AlO) 



(All) 
(A12) 



Here the distribution function n^ is considered as a matrix diagonal in momentum space, while the matrix $ has the 
matrix elements 



$fei-fc2 =B(fci -fcs). 



2tt 



:h,^ 



L{ki — fc2 
Using the fermionic commutation relations satisfied by Cfc and the fact that all Ck annihilate \N), one now finds 



(A13) 



exp 



E/5«^. 



,?>0 



|iV) = exp 



E«W^fc,; 



.kik2 



\N) 



W^ I (isexp[sf/]l/exp[st/'^] 



(A14) 



(A15) 



Note that the operator $ preserves the subspace k > ^N while the operator $'^ preserves the subspace k < ^N. 
On the other hand, on the first of these subspaces 1 — n^ — 1, while the other subspace is nullified by 1 — n^ . It 
follows that 



e^U = e'^(l-n'') + e-^'' n^^ 



Using this relation, the expression for the matrix W can be simplified to the form 

1 



W 



The matrix element (IA6I) now reads 



,-*' 



(1 - n^) e*" - e* (1 - n^) 



JV^g-* 



(A16) 



(A17) 



{l3*,N\p\N,(3)^j{N\exp 



Yck,Ck2W+^ 

.kik2 



exp 



z2 ^fei'^fe2^fci* 

kik2 



exp 



.ki k2 



\N) 



= l{N\e^^e^e^\N) , (A18) 



where 



Ek^M = e(fci) (l - 2n^(/ei)) 4^^^ 7 7 = exp 



J2eik){n^{k)~n'ik)) 



(A19) 



For the evaluation of the average (lAlSp one can employ the standard technique of the fermionic coherent states. Let 
us define the set of the eigenstates of annihilation operators Ck with the Grassmann eigenvalues Xk (x stands for the 
full set of Xk with k = —00 . . . 00): 



Cfclx) =Xfc|x)- 



(A20) 
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The coherent states \x) form an overcomplete basis with the resolution of identity 



1 = / Dx*Dxe'x.p 



k 



\x){x*\ 



where 



Dx*Dx^l[dxldxk 



(A21) 



(A22) 



Using the resolution of identity (|A2ip and noting that the state |iV) is itself a coherent state |x = 0), we get 



{(3* , N\p\N , /?) = I / Dx'*Dx'Dx'*Dx' exp 



-J2i^k'xl + xl*xl) 



X ix* ^0\e^' \x'){x'*\e'^\x'){x'*\e'^ \x = 0) . (A23) 



Finally, since 



{X*\e \x) =cxp 



E 



xlXk+ Yl XfeiXL^fci. 

k ki ,k2 



(X*|e^lx) =exp 



/ . XkiXk2 [^ Ikif, 

ki,k2 



we are left with a simple Gaussian integral over Grassman variables. The result of the integration reads 



{l3*,N\p\N,l3} = l 



det 



1 -2e-'^We-'^ 



1/2 



2W+ 1 

Using the explicit form of the matrices W and E, one can reduce the expression above to 



{p*,N\p\N,P) -^det e-^e* (1 - nw) e-^e-^ + e-* TiAre" 



Finally, working out the factor in front of the determinant in (IA28P , we find 

(/3*,iV|p|Ar, /3) = det [l - e*njve"*(l - n) - e-*^(l - nA.)e*^n 



(A24) 

(A25) 
(A26) 

(A27) 

(A28) 

(A29) 



Equation (jA29|) gives an explicit expression for the density matrix of non-equilibrium interacting fermions in the 
bosonic basis. It has a form of a one-dimensional functional determinant of Fredholm type. 

Let us verify that Eq. (IA29[) reduces to the Boltzmann-Gibbs form at equilibrium, which in present notations 
corresponds to e(fc) = kVpjT . Indeed, we know that for e(fc) = Eq. (|A28p should give just the norm of the coherent 
|7V, /3) (up to the normalization factor Z) 



(/3*,iV|p|7V,/3)|^^^^^^ 



■det 



e (1 — njss) e + e njqe 



<i)+ $+ 



: exp 



.g>0 



(A30) 



On the other hand for e(fc) == kVp/T we can write 



— <!>+ $+ 



(^*,A^|peg|A^,^) = ^det e''e*e'(l-njv)e~'e-*e' + e-* njye 



= ^det 



e (1 — Tijv) e + e riNC 



where 



*fei-fe2 



[^"^*^l....= 0(^1-^2 



27r 



Hh - fc2) 



/3fe,_fe,e-^-('=-^-=)/^, 



(A31) 
(A32) 
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We see that the matrix elements of the equiUbrium density matrix can be obtained from the norm of the coherent 
states just by the replacement 



/3,^/3,e-^-'/^. 



with the result 



(/3*,Ar|p^,|iV,/3) = -exp 



.i3>0 



m 



In the operator form, Eq. (jA34|) reads 



Peq = -^ cxp 



^E'^^J^. 



g>0 



(A33) 



(A34) 



(A35) 



which is the expected equilibrium result. 

Away from equilibrium, the bosonic density matrix (|A29p is much more complicated. Nevertheless, it can be used 
to derive the expressions (|T7lH51H^ for the prototypical average (H5)l required for evaluation of correlation functions. 
To do this, one represents the determinant (jA29[) as an a Gaussian integral over auxiliary Grassmann variables and 
performs an expansion in powers of e*. At every order, evaluation of the correlator of the bosonic exponents (j43p is 
then given by a Gaussian integral over bosons and fermions. The crucial point is some specific cancellations between 
the fermionic and bosonic integrations, which finally allows us to evaluate all orders of the expansion analytically and 
resum the whole series. This brute-force derivation requires, however, rather cumbersome combinatorics, and we do 
not present it here. An alternative, considerably simpler, derivation involves refermionization, as explained in the 
main text and Appendix IB] 

Appendix B: Averaging exponential of bosonic fields in a non-equilibrium Luttinger liquid via 

refermionization 



In this Appendix we evaluate the prototypical correlation function (|43p where the average is performed with the 
density matrix (|4ip. Since the right-moving in- fermions are completely independent from the left- moving ones (see 
Eq. (j41l) ) we ignore the latter for a while and include the left-movers in the final formulas. In our notations we also 
suppress the index "in" since all the operators (bosons and fermions) we will be dealing with in this Appendix are 
in-operators and the suppression should not cause any confusion. Thus, we need to evaluate the average 



Z[g\g] 



exp 



E3(«)^' 

L?>0 



exp 



q>Q 



with the density matrix 



P = ^ exp 



-Y,eik){ata,-n'ik)) 



(Bl) 



(B2) 



We can now apply the well-known expression for the trace over the fermionic Hilbert space of a product of exponentially 
quadratic operators^: 



tre 



ffi 



,ff„ 



det(l 



ah„' 



where 



H.. = Y: h^^'^al 



a-ko 



(B3) 



(B4) 



ki,k2 



The right-hand side of (jB3p is a determinant of an operator acting in the single-particle Hilbert space. Applying (jB3 
to the average in question, we get 



Z[g*,g] = det ( 1 - n{k) + e-'^'^'''>n{k) 
P2^ 



5{x) 



q>0 



{9q( 



,iqx 



*p-l<]X\ 



9qe 



(B5) 
(B6) 



22 

The determinant in (jB5[) is in fact not well defined, for the following reason. Strictly speaking, Eq. (jB3|) assumes 
that in the fermionic Hilbert space there is a state nullified by all the operators Ofc . This is not the case in the present 
situation. We can resolve this difficulty by making the transformation from the particle operators a^, a^ to particle 
and hole operators: 

Cfe=afc(l-n°(fc))+a+n°(fc), (B7) 

4 =at{l-n"{k))+akn"{k). (B8) 

Then we can evaluate Z[g*,g] with the technique of fermionic coherent states. The calculation is very similar to that 
of Appendix lAl 

We can make a short-cut, however, if we notice that Z[g* , g] is an average of an operator normal-ordered in bosons. 
Thus, Z[g*,g] = 1 at zero temperature. This suggests that the regularized version of Eq. (|B5[) should be 



Z[g\g] ^det 



(l - n°{k) + e^'''(^)n°(fc)) ^ (l - n{k) + e~'^^''^n{k)') 



(B9) 



A direct calculation along the lines of Appendix 1X1 indeed confirms this result. 

Incorporating now the left electrons into our consideration, we come to Eqs. ([47l [48l |49)) . 
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